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Abstract. We introduce and discuss an effective model of a self- gravitating system 
whose equilibrium thermodynamics can be solved in both the microcanonical and the 
canonical ensemble, up to a maximization with respect to a single variable. Such 
"j^ , a model can be derived from a model of self-gravitating particles confined on a 

' ring, referred to as the self-gravitating ring (SGR) model, allowing a quantitative 

comparison between the thermodynamics of the two models. Despite the rather crude 
approximations involved in its derivation, the effective model compares quite well 
with the SGR model. Moreover, we discuss the relation between the effective model 
, presented here and another model introduced by Thirring forty years ago. The two 

^ ' models are very similar and can be considered as examples of a class of minimal models 

^ ■ of self-gravitating systems. 

> ■ PACS numbers: 05.70.Fh; 95.10.Ce 

' 1. Introduction 

2 . Systems of classical particles mutually interacting via gravitational forces can model the 
O . behavior of many objects in the universe (globular clusters, elliptical galaxies, clusters 
of galaxies) as long as other interactions are negligible compared to the gravitational 
ones [1]. When the number of particles is large the direct numerical simulation of 
such systems is a heavy task [2] and it would be reasonable and useful to approach them 
via equilibrium statistical mechanics. However, self-gravitating systems do not have 
a "true" thermal equilibrium for two main reasons [3]: (i) the gravitational potential 
is singular for vanishing distance between two particles, making (at least part of) the 
system collapse in states with infinite density and (ii) particles that do not collapse tend 
to escape the system (evaporation). From a physical point of view the first problem can 
be easily solved. No real system exists where the only non-negligible interaction is 
classical gravity at all length scales: either the interacting "particles" are macroscopic 
bodies like stars or galaxies, or quantum effects must be taken into account below a 
certain length scale. In both length scale exists (the size of the bodies or 

the scale where quantum effects set in) below which the potential has no longer the 
classical gravitational form. If one is not interested in small-scale details, the potential 
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can be regularized by replacing it with a softened one at short distances (U E] or by 
directly considering self-gravitating fermions j5]. To solve the second problem is less 
straightforward and one is forced to (somehow artificially) confine the particles in a 
finite volume. However, on physical grounds such an approximation is reasonable since in 
many cases the evaporation rate is slow compared to the other time scales in the system 
[6] . A regularized and confined self-gravitating system has a thermal equilibrium in both 
the canonical and the microcanonical ensemble [7J. Such a system can thus be studied 
within the framework of equilibrium statistical mechanics, although its behaviour in 
the two ensembles is very different: it can be considered a prototype of systems with 
ensemble inequivalence [8], showing e.g. a core-halo phase with negative specific heat in 
the microcanonical ensemble which is replaced by a discontinuous phase transition from 
a clustered to a gas phase in the canonical ensemble |9]. 

Canonical and microcanonical MonteCarlo simulations of a full self-gravitating 
system in three spatial dimensions are heavy [9J , although they may be a little easier than 
the direct integration of the equations of motion. This suggests to look for simplified 
models which may be easier to study; one main simplification comes from considering 
models which are effectively one-dimensional. Many such models have been introduced 
in the last decades, ranging from the sheet model [10] and the shell model [11] to the 
self-gravitating ring (SGR) model [12]. The latter is particularly interesting because 
the interaction among the particles is given by the full three-dimensional gravitational 
potential (regularized at small distances), while the particles are confined on a ring. This 
yields a behaviour that is qualitatively very close to that found in three-dimensional 
systems, although allowing a much simpler study. In the limit of an infinite number of 
particles, the model can be studied in the mean-field approximation with a very efficient 
numerical technique [13], showing that in the microcanonical ensemble there is a phase 
transition separating a homogeneous high-energy phase from a clustered phase. An 
independent analytical argument supporting the existence of such a transition has been 
given in [I^ . 

The aim of the present paper is to introduce and discuss an effective model which 
approximates the SGR model and is exactly solvable (up to a maximization with 
respect to a single variable which has to be performed numerically). Although the 
approximations involved are rather crude, the behaviour of this effective model is very 
close to that of the SGR model, not only at a qualitative level but also at a semi- 
quantitative one: this suggests that the approximations made do capture most of the 
important physics. 

The paper is organized as follows. After recalling the main features of the SGR 
model we introduce and discuss our effective model, and we present its solution in 
both the microcanonical and the canonical ensembles; the details of the calculations 
are reported in two appendices. Then, we compare our model with one introduced 
by Thirring forty years ago [15j: the two models are indeed very similar, although the 
latter was not aimed at approximating any particular explicit model. We end with some 
remarks and a discussion of some open issues. 
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2. Model 

We now want to introduce and discuss our model. Before doing that, we briefly recall 
the main features of the self-gravitating ring (SGR) model. 

2.1. Self- gravitating ring model 

The SGR model describes identical classical particles of mass m constrained to move 
on a ring of radius R and mutually interacting via a regularized gravitational potential. 
Its Hamiltonian is [12] 



^SGR 2^^2 ^ 2 2^^ v^i?v/l-cos(^.-^,) + a ' ^ ^ 

where di G (— tt, tt] is the angular coordinate of the z-th particle, 

P. ^ (2) 

is its angular momentum, R is the radius of the circle and G is Newton's gravitational 
constant. The constant a > provides the short-distance regularization: when 
di — dj'^l the interaction between the i-th and the j-th particle is effectively harmonic 
on a length scale of the order of 

d„ = R\/2a . (3) 

The SGR model can thus describe a self-gravitating system when da ^ R, i.e., a <^1. 
On the other hand, in the opposite limit a — )■ oo the SGR model becomes equivalent 
to the ferromagnetic hamiltonian mean-fleld (HMF) model [16], that is a model of fully 
coupled planar classical spins with ferromagnetic interactions. 

In the microcanonical ensemble, the fundamental quantity to compute is the density 
of states. For the Hamiltonian (JT]) it can be written as 

UJn{E) = — dPi---dPM ddi---ddN5{'^SGK-E) . (4) 

' • J —oo J — TT 

Deflning a characteristic time scale r as 

R3 

(5) 



GNm 

and dimensionless momenta as 
dd 

Pi = —, t = l,...,N (6) 
one can write a dimensionless Hamiltonian "H as 

TV N 

'H = -yp^-^-^y ^= ^ (7) 



1=1 * *ii=i * V * J 

Since 

'H = — ^'HsGR ; (8) 
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the density of states (jl]) can be written as 



ujn{E) = — Wm~ dpi---dpN ■ ■ ■ 5 

T -1-^ ■ J — oo J — IT 



n-E 



,(9) 



~ 2 

where E = -^^^E, and using the properties of the 5 function it becomes 



^7V-1^2(JV-l) 



^^(^) = —p^r2Y^^N{E) (10) 

where a) at is the dimensionless density of states 

Cj^{E) = I dpi---dpN I d^i--- d^N S(n-E) . (11) 

J — oo J —TT 

Before going on, let us note that the entropy derived from the full density of states 
is extensive in the thermodynamic limit A^— J-oo,/?— j-oo, ^— )■ const.; moreover, 
the adimensional Hamiltonian ([7]) is also extensive, due to the rescaling of the 
coupling between the particles. The Kac prescription for making extensive a long- 
range interaction [S] is here obtained via a suitable choice of the adimensionalization 
procedure. 

From now on we will consider only the dimensionless Hamiltonian ([7]) and density 
of states ( fTTl) . For ease of notation, we will drop the tilda's and simply denote them by 
Ti and un{E), where also E is adimensional. 



2.2. Effective model 

We now want to approximate the Hamiltonian of the SGR model, or more precisely the 
adimensional Hamiltonian ([7]), in order to make it soluble, i.e., such that the density of 
states (ITT]) can be computed. 

Numerical simulations of the dynamics of the SGR model reported in [12] have 
shown that, at a given energy, particles can be roughly divided into three classes 
according to their dynamical behaviour: cluster particles, gas particles, and halo 
particles. Cluster particles are tightly bound in a cluster and never get far from it; 
gas particles move almost freely around the ring; halo particles have a complicated 
dynamics that is somehow intermediate between the other two. The relative population 
of particles in the three classes depends on the energy (or temperature): at low energy 
almost all the particles are cluster particles, while at high energy all the particles are 
gas particles. 

The strategy we are going to implement in order to define an effective model is to 
consider only the first two classes of particles (cluster and gas) and to assume that each 
particle belongs to one of the two classes. This allows a simplification of the potential 
energy which makes the model soluble. 

Let us then assume that Ng particles, with 1 < Ng < N, are gas particles. Then, 
writing the Hamiltonian as 

1 ^ 

n = -J2p' + v{A,...,^N), (12) 

1=1 
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we can split the potential energy V into three parts: 

+ V,^t{A,...,^N) , (13) 



where 



where 



1 ' 

i,jr = l 
1 ^ 

Cluster (^jV,+l, . . . , ^jv) = - ^^r^ (15) 

i,j=Afg+l 

i=l j = Ng + l 

= „ ^ ^ ■ (17) 
Vl — cosx + a 

Up to this point we have only rewritten the potential energy in a different form. However, 
this form naturally allows to introduce the approximations which make the model 
soluble. Let us now discuss the approximations. 

{i) Since the Ng gas particles are essentially free particles, as far as their interaction 
V^as is concerned we consider them as uniformly distributed on the circle. The 
interaction energy f|T^ between these particles is then a constant: 

Kas = -7 ^ (18) 

^ '2NV2 

where 

,^^/\,W..^-|=,c(j|^) (19) 

and /C(x) is the complete elliptic integral of the first kind. 

[a) We consider the N — Ng remaining particles as confined in a cluster. We assume 
the cluster is tight, i.e., the particles are all close to each other: 

- 1?,- < 1 = Ng + 1,...,N . (20) 

We can thus expand the interaction energy (1151) among these particles up to the 
harmonic order, and write 

Cluster = + — E (^^ - ^^)' 5 (21) 



such an approximation is reliable if an assumption stronger than (1201) holds, i.e. 



\2 



<1 Wi,j = Ng + l,...,N . (22) 



a ' 
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Moreover, since the particles in the cluster do not "feel" the topology of the 
circle, we assume that 

-oo <^i< +00, yi = Ng + 1, . . . , N ; (23) 

this will allow the analytical computation of the configurational integrals in the 
density of states. 

{Hi) As far as the interaction f|T6|) between cluster and gas particles is concerned, we 
note that as long as the assumption ( l20l) holds, the typical distance between a gas 
particle and a cluster particel is much larger than typical interparticle distances in 
the cluster, so that we may assume that all the cluster particles are in the same 
location, i.e., = 0. Being the gas particles uniformly distributed on the circle, 
this yields a constant for l^nt, i-e., 

where 7 is given by f[T^ . 
The Hamiltonian of our effective model is then 

1 ^ 

■^eff = ^5Zp-+KfI , (25) 



where 



2 

i=l 



N 

Ks = -Vo{N,Ng,a) + ^ J2 i^^-^J)\ (26) 

i,j = Ng + l 



and where we have set 

\2 



f/V - /V V 



9 



N\/2 2NV2, 



(27) 



and 



^ 2iV(2a)3/2 ■ 



3. Microcanonical and canonical thermodynamics 

Let us now discuss the solution of the effective model in the microcanonical and canonical 
ensembles. In the limit — > 00, at fixed Ng the model is exactly solvable in both 
ensembles. However, Ng is not a priori assigned and must be fixed in a self-consistent 
way. The simplest way to do so is to take into account all the possible values of A^^; 
as we shall show in the following, in the limit — )■ 00 the model is still solvable up 
to a maximization (resp. minimization) in a single variable which must be performed 
numerically, whose physical meaning is just to determine the value of Ng that maximizes 
the entropy (resp. minimizes the free energy). 
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3.1. Microcanonical ensemble 



To solve the model in the microcanonical ensemble we need to calculate the entropy 
density 



s{e) = lim —logUNie) , 

N-¥oo iV 



(29) 



where we have introduced the energy density £ = and a; at is the density of states f lTT]) 
where "H is replaced by Ties'- 

/oo 
dpi--- dpN 
-oo 



d^i--- dd 



No 



dd 



d^N 5 CHes - Ne) 



N 



E 



N„=0 9 



Ng\{N-Ng)\J_ 



dpi ■ ■ ■ dpj\f / d-di ■ ■ ■ 



No 



(30) 



X 



^ N N 



d'&N.+i ■ ■ ■ d'dN S ^ ^ 

i=l ~ i,j = Ng + l 

In the above expression we have summed over all the possible choices of Ng, properly 
counted by the degeneracy factor ) . 

The calculation of the above integral is straightforward, albeit a bit involved, and 
can be performed following a procedure similar to that used in pTf [T8]. The details are 



reported in [Appendix A[ It turns out that the entropy density in the thermodynamic 
limit is given by 



s{e) = sup s{e, Ug) . 

"s6[0,n'^='>=(£)] 



N, 



/ and n^^^l 



(31) 



e) is the 



where we have introduced the fraction of gas particles ng 
maximum fraction of gas particles allowed at a given energy density e, given by Eq. 
f lA.lSp . The explicit expression of s{e,ng) is 

47r(2a)3/2 



Tin 



s[e,ng) 



+ 



2 

2 - n. 



log 



(l-%)(2-%). 
^ [1 + log a{ng, a, e)] + Ug log(2^) 

- Hg log rig - (1 - Ug) lOg ( 1 - U g) , 



2n^/2 
2 - n„ 



(32) 



with 



a{ng, a, e) 



^ ng{2-ng)+^ 



2^2 



2x/2a 



(33) 



1 



one can check that aijig, a, £:) > if ng G [0, n^^^{e)] and e > emm, where Emin ■ 
is the absolute minimum of the potential energy per degree of freedom. 

As anticipated, the solution of the effective model in the microcanonical ensemble 
amounts to finding the value n^{s) of ng realizing the extremum in (IHT]) . This can be 
easily done numerically, since the explicit form f p2|) of s{e, ng) is available. 




Figure 1. Effective model in the microcanonical ensemble, a = 10^^. (top) Fraction 
of gas particles rr^{s) (red line); (bottom) temperature T(e) computed for the effective 
model (red line) and for the SGR model (blue symbols). 



3.1.1. Results for the thermodynamic quantities In the following we report the results 
for the fraction of gas particles n^{s) and for the caloric curve, i.e., the temperature 
T{e) = (^) , as a function of e. We also compare the latter quantity with that 
computed for the SGR model via the numerical method introduced in [13]; for n^{e) 
such a comparison is impossible, because no such quantity is easily defined for the 
SGR model. In Fig. [T] we report n^{s) and T{e) computed for a softening parameter 
a = 10~^, as well as a comparison with T{e) for the SGR model; in Fig. [2] we report the 
same quantities for a = 3 x 10~^. The agreement with the SGR model is reasonably 
good already at a = 10~^ and becomes very good at a = 3 x 10~^. In both cases we find 
a phase transition from a homogeneous phase (characterized by ~ 1) to a clustered 
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20 



T 




20 



Figure 2. As in Fig. [T]for a = 3 x 10" 



phase while lowering e below a critical value Sc, the critical energy is Ec — —0.46 for 
a = 10~^ and Ec — —0.8 for a = 3 x 10~^. These values should be compared with those 
found for the SGR model, i.e., Ec — —0.32 for a = 10~^ and Ec — —0.5 for a = 3 x 10~^. 
The agreement is good, especially for the lower value of a. In the case a = 10~^ 
the phase transition is continuous, while it is discontinuous (the temperature T jumps 
between two different values at ^c) at a = 3 x 10~^. We find indeed a microcanonical 
tricritical point which is located at a ~ 5 x 10"'^; in the case of the SGR model this 
point is located at a ~ 10~^, so that also in this respect the two models are very similar. 
We stress that the above results come from a numerical maximization so that we have 
no rigorous proof of the existence of true singularities in the microcanonical ensemble 
for the effective model. However, the same holds for the SGR model too, as long as the 
softening parameter a is finite. 
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In the SGR model, the hirfi-energy phase {e > Ec) is a perfect gas phase, as shown 
by the numerical calculation^jl reported in [13]. We should thus expect = 1 for 
6 > Ec- This does not happen, although is very close to 1 (W^ ~ 0.93 for e > Ec when 
a = 10~^ and ~ 0.995 for e > Ec when a = 3 x 10~^). This is due to the presence of 
the degeneracy term (^) in the density of states, which makes the entropy vanish for 
rig = 1 and = 0, so that the extremum of the entropy fl32l) can never be realized in 
the boundaries of the domain of Ug when the energy density is strictly larger than its 
absolute minimum Emm- For the latter value of e, however, we have = 0, and there 
is a whole region of energy densities where ~ 0. This can already be seen in the top 
panels of Figs. [T] and [H however, it is more evident if we look at the first and second 
derivatives of %(£), reported in Fig. [3] in the case a = 3 x 10~^. The peak in the second 
derivative of n^{E) can be effectively taken as the upper limit of a "highly clustered 
phase" (which is not a "true" phase because there is no singularity). The existence 
of such a region is a very interesting feature of this model: it reminds what has been 
observed in simulations of confined and regularized three-dimensional self-gravitating 
systems [H |9], where one finds a low-energy clustered phase, an intermediate-energy 
"core-halo" phase (which would be mimicked by the coexistence of gas and cluster 
particles in the effective model) and a high-energy perfect gas phase. We shall come 
back to this point in the two final sections of the paper. 

An unphysical feature of the model is that when e is very large the number of 
gas particles starts to decrease and eventually ri^ — )■ for e — oo. This happens for 
any value of a and is due to the fact that the coordinates of the cluster particles are 
allowed to take values in M instead of in (— 7r,7r], so that at very high energy it is 
always convenient to make a "loose" cluster whose effective size is larger than the circle. 
However, this happens at values of e which get larger and larger as a gets smaller, so 
that one can safely ignore this fact in practice: in Fig. [T]one can see a first hint of this 
phenomenon, which is instead invisible in the energy range of Fig. O 

A feature of the model that does not compare well with the SGR model is that 
for any value of a there is a region of nonconcave entropy: there is always ensemble 
inequivalence, so that, as we shall see in the following, there is no canonical tricritical 
point and the phase transition in the canonical ensemble is always discontinuous. In the 
SGR model, on the contrary, the ensemble inequivalence is present only for a < ckct 
with ctcT — 0.1. On the other hand, the approximations made to derive the effective 
model are only reasonable for small values of a, so that this is not a big surprise. 



3.2. Canonical ensemble 

Let us now discuss the solution of the effective model in the canonical ensemble. We have 
to compute the partition function: the calculation is even more straightforward than in 



the microcanonical case, so that we report the details of the calculation in Appendix B 



I In the limit a — > oo, when the SGR model becomes equivalent to the ferromagnetic HMF model, this 
can be shown analytically; however, for large values of a our effective model is not a good approximation. 
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Figure 3. Effective model in the microcanonical ensemble, a = 3 x 10 ^ . (top) First 
derivative of the fraction of gas particles, (bottom) Second derivative -j^- 



It turns out that in the thermodynamic hmit N ^ oo the free energy density is given 
by 

/(/?)= inf f(/3,n,) (34) 

nge [0,1] 

where /3 = T^^ (we set the Boltzmann cosntant to unity) and 



log(l -Ug)-^ log(27r) - — log 



2n 

J 



ng'yirig - 2) (1 - ^-rig 



log 



m - ng) 



(35) 



2^2 2V2a 2/3 

Again, the solution amounts to finding the value n^{P) of Ug realizing the extremum 
in (^^. This can be easily done numerically, since the explicit form (^5^ of f{f3,ng) is 
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Figure 4. Effective model, a — 10^^. Comparison between the caloric curve 
T(ii) computed in the canonical ensemble (blue line) and that computed in the 
microcanonical ensemble (red line; already reported in the bottom panel of Fig. [1} . 
On the horizontal axis there is u in the canonical case and s in the microcanonical 
case. 



3.2.1. Results for the thermodynamic quantities In Fig. IHwe report the results for the 
energy density u = ^^^4P-, i.e., 



for a = 10~^ and we compare them with the microcanonical results. In order to make 
the comparison easier, we do not plot m(/3) or m(T); we rather report T as a function of 
u, i.e., the canonical caloric curve. 

As already noted, there is ensemble inequivalence for any value of the softening 
parameter a. Moreover, the region of inequivalence is larger than in the SGR model: 
for instance, when a = 10~^ the inequivalence region is bounded below by an energy 
density miow — —3.15, while in the SGR model the bound is at miow — —1-98. 

4. Thirring model 

We now discuss the relation between our model and a model introduced by W. Thirring 
forty years ago [I5] as a toy model of a system^ with negative specific heat. The latter 
model is solvable up to a maximization in a single variable, as our effective model. With 
a suitable choice of the parameters the microcanonical thermodynamics of the Thirring 
model is qualitatively similar to that of self-gravitating systems, so that we expect it to 

§ Thirring called it "A somewhat artificial version of a star." 



known. 




(36) 
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be close to that of the effective model introduced above. It will turn out that the two 
models are even more closely related than one would a priori suspect. 

To define the Thirring model, consider classical particles confined in a three- 
dimensional volume V and a volume Vq G V, which will play the role of the cluster. 
The Hamiltoniar[ii] is 

TV N 

^=25:P^' + ^E^(r^'r.) (37) 
1=1 «j=i 

where the potential energy is highly nonlocal, i.e., v is such that particles i and j interact 
with constant energy if they are both in Vq and are free otherwise: 

v{ri,rj) = -2iyxvo{ri)xvo{rj) , (38) 
where ly > (the force is attractive) and xvo is the characteristic function of the subset 
Vq, i.e., XVoi^) = 1 if a; G Vo and XVoi^) = if x ^ Vq. Following [15], we write the 
volume outside Vq as 

V-Vo = e^Vo ; (39) 

when e^ ^ 1, i.e., the volume outside Vq is much larger than Vq, in the microcanonical 
ensemble the model is solvable in the thermodynamic limit — )■ oo up to a maximization 
on a single variable, the number Nq of particles in Vq- Such a maximization comes from 
the evaluation of the density of states containing all the possible values of A^o via a 
saddle-point approximation, similarly to what we have done for the effective model. 
Details on the solution can be found in [15]. The entropy density can be written 



s{e) = sup 5(5,77-0) (40) 
noeW",i] 

where e = ^ is the energy density, no = ^ is the fraction of particles in Vq, n™™ is the 
minimum value of no allowed at a given e and 



s{£,no) = - log I — ) + log 



(£ + z/n2)3/2e^(i-"o)+i 
n"o(l -no)i-"o 



(41) 



up to an irrelevant additive constant. In the following we shall denote by no the value 
of no that maximizes the entropy (HTl) . 

Figure [5] shows an example of the thermodynamics of the Thirring model; we report 

I — no (to make the comparison with the effective model easier) and T as a function of 
e. Apart from the obvious differences in the slopes of the low- and high-energy parts of 
the T{e) curve, due to the different dimensionality and to the different regularization 
(free particles in a box vs. harmonic forces), the caloric curve T{e) of the Thirring model 
is remarkably similar to that of the effective model. This overall similarity is not a big 
surprise: it is a confirmation that the qualitative behavior of a confined, regularized 
self-gravitating system can be captured by a "cluster + gas" model regardless of the 
details on the regularization. These models are in a way the minimal models of the 
equilibrium statistical mechanics of self-gravitating systems. What is more interesting 

II Note that, at variance with Ref. |15j . we used the Kac prescription to make the Hamiltonian extensive. 
This will allow us to study the model in the "conventional" thermodynamic limit. 
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Figure 5. Thirring model in tlie microcanonical ensemble, with j/ = 1 and F — 6. 
(top) Fraction of particles outside Vq, 1 — rTo{e); (bottom) temperature T(e). 



is that this similarity extends to the low-energy behaviour, where one would expect 
a stronger dependence on the details of the regularization. Despite the fact that the 
small- distance regularization is completely different in the two models, Fig. |6] shows that 
the behaviour of the first and second derivatives of 1 — fiois) in the Thirring model is 
strikingly similar to that of n^{s) in the effective model (Fig. [3]), with a strong peak in 
the second derivative marking the upper limit of the strongly clustered regime, where 
riQ ~ 1. As in the effective model, due to degeneracy factor in the density of states there 
are no phases with tTq = 1 ot tTq = 0; a.t high energies, in the gas phase, no is very small 
but always strictly positiv^. 

% Due to the fact that the cluster size is fixed the Thirring model does not show the unphysical decrease 
of the fraction of particles outside the cluster at very high energies. 
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Figure 6. Thirring model in tlie microcanonical ensemble; numerical values of the 
parameters as in Fig. \E\ (top) First derivative of the fraction of particles outside Vq , 
(bottom) Second derivative 



At variance with the Thirring model, the effective model we have introduced 
above is derived from a "microscopic" model, the SGR model. This feature allows 
a detailed comparison between the thermodynamics of the effective model and that of 
the "microscopic" one, allowing to test the effect of the approximations. Apart from this 
feature, our effective model is very similar to the Thirring model, and can be considered 
as a way to derive "Thirring-like" (i.e., "cluster + gas") toy models from models with 
true gravitational interactions. 
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5. Conclusions and outlook 

We have discussed how to derive an effective model of a self-gravitating system starting 
from the SGR model introduced and studied in Refs. [121 [13]. Such an effective model 
is solvable up to a numerical maximization in a single variable (which plays the role of 
an order parameter) in both the microcanonical and the canonical ensembles. Following 
a suggestion coming from the study of the dynamics of the SGR model [121 US], the 
effective model assumes that particles can be split into two families: cluster particles, 
all of which are interacting with each other by harmonic forces, and gas particles, which 
feel only a constant potential due to both the cluster particles and the mean field of 
the other gas particles. The fraction of gas particles, Ug, is the "order parameter" to be 
determined self-consistently. 

Despite the rather crude approximations used to derive the effective model, the 
results for the thermodynamic quantities are quite close to those found for the SGR 
model using the numerical method developed in [I3], especially for small values of the 
softening parameter a, where the behaviour of the two systems is closer to that of 
an "ideal" self-gravitating system. Although one could expect an agreement at small 
energies, the agreement is definitely good also at energies up to and above the transition 
between the phase with negative specific heat and the homogeneous phase, which is well 
reproduced by our results. It is interesting that such a simple toy model as ours is able 
to reproduce the thermodynamics of the SGR model even at quantitative level. 

A qualitative disagreement appears at larger values of the softening parameter, 
where in the SGR model the entropy becomes concave and there is no longer ensemble 
inequivalence: in our effective model the entropy appears to be always nonconcave, for 
any value of a; however, the approximations made are no longer justified when a is not 
small. 

Clearly, there is room for improvements: the caloric curve T{e) is well reproduced 
when a is small, but its shape in the low-energy part — see e.g. Fig. |2] — shows that there 
are important anharmonic effects which should be taken into account. 

We have also compared our model with another model introduced by Thirring in 
1970 [15]. In that model one has a cluster of particles confined in a finite region of space, 
all interacting with each other via a constant attractive potential. This region of space 
is enclosed in larger volume where particles can escape becoming free particles. When 
the cluster extension is small, the thermodynamics of the two models is very similar, 
also in the low energy region dominated by the regularization of the potential, although 
the nature of the cluster is different in the two models. This shows, in our opinion, 
that "cluster + gas" toy models like these are good candidates as minimal models of 
self-gravitating systems. 

An interesting feature of our model (and of the Thirring model too, although this 
had not been noticed before, to the best of our knowledge) is the presence of a low-energy 
region where the fraction of gas particles rig is very small and stays very small up to 
a certain energy where it starts rising rapidly. There is a mathematical reason why rig 
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can not be exactly zero in a finite region of energy, i.e., tlie degeneracy factor associated 
with counting tlie number of ways of constructing a state with a given fraction of gas 
particles. This counting assumes that all the gas particles are independent, which is 
clearly an oversimplification. One could wonder whether a more refined counting may 
imply the presence of a singularity at low energy, bounding a phase with Ug = 0, whose 
existence for the SGR model has been conjectured in [H]. A hint in this direction comes 
from assuming that the gas particles are all correlated like bosons (which is wrong, but 
is somehow the opposite situation to that considered here), so that states obtained from 
each other by interchanging two gas particles should not be considered as distinct and 
no degeneracy factor would be present in the density of states. Such a calculation has 
been performed in [19] and yields a sharp transition at a finite energy density and a 
phase with Ug = 0. The drawback is that the agreement with the thermodynamics of the 
SGR model is definitely worse. It is tempting to think that maybe a weaker transition 
occurs, and that it can be described by a "proper" counting of the degeneracy of the 
gas particle states. 
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Appendix A. Density of states and entropy 



To solve the model in the microcanonical ensemble we need to calculate the density of 
states (1301) . i.e.. 



N 



Ng\{N-Ng)\ 



dpi ■ ■ ■ dpN I d'di ■ ■ ■ dd 



Na 



X 6 



N 



N 



i=l iJ=Ng + l 

To compute the above integral we follow a procedure similar to that used in pTl [T8] 
First, we expand the square in the last sum, obtaining 



d-dNg+i ■ ■ ■ d^&N 
(A.l) 



UN 



1 ^ 



iV! 



dpi ■ ■ ■ dpN / di!}i ■ ■ ■ dd 



dd 



N„+l 



■dd 



N 



X 6 



N 



i=Ng + l 



A = Ng + l 



Vq + E 

1^ 



(A.2) 



and then we search for a coordinate transformation diagonalizing the coupling between 
the '(9's. As to this point, we note that 

2 



N 



N 



(N-Ng) Y ^^-\ E ^0 =i^Ng + l,---,^N)M^Ng + l,...,^Nf 



i=Ng+l 



. i=Ng+l 
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(A.4) 



where the symmetric {N — Ng) x [N — Ng) matrix A reads as 

/I ■■■ 1\ 

A = - ; ••. ; +{N-Ng)iN-N, 
\i ... i) 

and 1(1 is the d x d identity matrix. The matrix A has eigenvalues Ai = and 
A2 = ■ ■ ■ = ^N-Ng = N — Ng, and can be diagonahzed by an orthogonal transformation 
that does not change the integration measure. Hence we can write 

N 



UN 



X 6 



Ng\{N-Ng)\ 



{2n 



dpi . . . dp 



N 



dd 



N„+l 



dd 



N 



^ i=l i=Ng+l 



(A.5) 



where we have also performed the Ng integrals over the circle. With the change of 
variables ijji = di^jN — Ng, i = Ng + 1, . . . , N — 1 we get 

N 



UN 



X 6 



Ng\{N-Ng)\ 



dpi . . . dp 



N 



N N-l 
^ 2=1 i=Ng+l 



Vq + E 



N 



The last integral in Eq. ( 1A.6I) stems from the zero mode due to the 0(2) invariance of 
the Hamiltonian; it is divergent but does not affect the thermodynamic quantities so 
that from now on we will ignore it. 

2(Vo + E) — Yl!i=iP'i — integrals over the V^'s give the volume of 

2{Vo + E)-ZIiP1 



When f 



the {N — Ng — 2)-dimensional sphere ^® ^ of radius i? = 77- 



on the other hand, when 77- 

' 2fj, 



2{Vo + E)-J:Lp\ 



2/x 



< the same integrals vanish. As 



to the integrations over the momenta, we note that the integrand depends only on 



p 



N 



(A,7) 



SO that we can write 



N 



X 



2^/^ / dpp''-' 



Ng\{N-Ng)\ 

2^' -p2n {N-Ng-2)/2 



p , N-Ng -1 



r(f) 



2/i 



e 



2E'-p^ 
2JI 



(A.8) 



where T{x) is the Euler gamma function, Q{x) is the Heaviside step function and 
E' = Vo + E. 
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We want to compute 

s{e) = lim ^logUNiNe) , (A.9) 

where we have introduced the energy density e = ^. Clearly, 

ojN{Ne) = Oif e<e^i^, (A. 10) 

where e^i^ is the absolute minimum of the potential energy per degree of freedom, 

(A.ll) 



2V2a 

so that the domain of the entropy density (IA.9P is e > emm- For large A^, computing 
the integral in Eq. f lA.Sp with the Laplace (or saddle-point) method [20] we get 

/\rmax('_\ 



2N/2 

UNiNe) = V (2 



N,\{N-N,)\ 



3/21 (7V-iV9-2)/2 



X [iV22(2a)3/2] 
2 — n 



exp <|^A^ 
1 — n 



^log 



r 



2 -nr. 



r(f) 

(A.12) 



log 



2 - 



where 



aK,a,.) = ^n,(2-n,) + ii-^ 



(A.13) 



and we have introduced the fraction of gas particles Ug = N^^^{e) is the maximum 
number of gas particles allowed at a given energy density e, so that the quantity 
a{ng, a, e) given by Eq. flA.lSP is positive in the domain e > Emm with < A^g < N^^^. 

Neglecting the sub-exponential terms and using the Stirling approximation, Eq. 
flA.12p can be written as 

dngexp[Ns{e,ng)] , 



wher^ 



uJN{Ne) 



n--(£) 







and 



1 -n„ 



1 - 



1 + 2£V2a 
^/2-'y^/a 



(A. 14) 



(A.15) 



s[e,ng) 



2 

2 - n, 



log 



47r(2a 



,3/2 



1 -%)(2 -%) 
-1 ^ [1 + log a{ng, a, e)] + Ug log(27r) 

- Ug log rig - (1 - Ug) 10g(l - U g) ] 



271^2 

2 - Tin 



hence, we can write the entropy (]A.9P as 
s{e) = sup s{e, Ug) . 

+ Eq. (jA.15[) holds for sufficiently small a, i.e., such that \/2 — "f^/a > 0. 



(A.16) 



(A.17) 
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We want to compute the partition function of the effective model. As in the 
microcanonical case, up to a factor that takes into account the correct dimensions and 
makes the free energy extensive, we can consider a dimensionless partition function: 

iv^ A! r 

= 1^77777 ^TTT / dpi--- dp^ / d'di--- di!)N, 



X 



A,!(A-A,)! 



d^M+i ■ ■ ■ d^N exp (-/3'Hefr) 



(B.l^ 



Using the expression (1251) for T-Les and performing the integrals over the momenta and 
the first A^ angles we can write 

m = CiY' E 1^7T7#^(2vr)-« exp { /3 



X 



B J ^ A '(A - A„)! 



d^Ng^i ■ ■ ■ d^N exp 



■^^ A,(2A-A,) ^ (A 




(B.2) 



where 



/3 



' (2a)3/2Ar ■ 
Using the Hubbard-Stratonovich formula 



dyexp {—ay^ + 6y) 



(B.3) 



(B.4) 



with a = 1 and b - 
gets 



2rj XlilAT +1 and performing the integration over the t^'s one 



27r 



N/2 N 



A! 



(27r)^« exp <; p 



Ng=0 9 



Ng\{N-Ng)\ 



X 



27r 



1 {N-NJ2) 



dy 



\ N,{2N-N,) ^ (A - 



2AV2^ 



(B.5) 



The integral in ( IB.SP is divergent (as in the microcanonical case, it comes from the zero 
mode associated to the rotational invariance of the potential) but does not affect the 
thermodynamic quantities and will be discarded from now on. 
The free energy density is 



/(/3) = -^ lim ^\ogZ{P). 

p Af-5-oo iV 



(B.6) 



In the limit A — )■ oo, the sum over A^ in (IB.SP is dominated by its largest term; using 
the Stirling approximation and introducing the fraction of gas particles = ^ we can 
thus write 



/(/?)= irif /(/3,n,) 

nse[0,l] 



(B.7) 
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where 




n,)--log(2x)--log ^- 



) 



1 - n„ , r27r(2a)3/2 

log — ^ — - — 

2/3 ^ L/3(l - n,) 
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